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Abstract 

The translational order in a three-dimensional lattice pinned by hopping defects is studied. 
It is suggested that the equilibrium state of the lattice crosses from a glass to a crystal on 
the rate of hopping. It is argued that such a transition may exist in flux lattices of high- 
temperature superconductors due to the hopping of oxygen. 



PACS. 74.60.Ge - Flux pinning, flux creep, and flux lattice dynamics. 
PACS. 05.20.-y - Statistical mechanics. 



It is known that below four dimensions the translational order in a lattice subject to a static 
random background (pinning) is destroyed by a however weak interaction with the background 
[[]], @, |3|, 3 • This fact has received much attention mostly in connection with pinning of flux lattices 
in superconductors, where translational correlations determine such important quantities as the 
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resistivity and the critical current. The problem has been re-examined lately on the issue of the 
explicit form of the correlation function at large distances [|, ^, [7], |S|, [T(J. Apart from that 
question, it is generally agreed that the correlation length of the translational order scales with the 
strength of pinning / as \/ f 2 K A ~ d ) ^ Thus, pinning, no matter how weak, leads to a glass 



state [11] characterized by short range translational correlations. This effect of the static disorder 
is very different from the effect of the thermal disorder, as one can see from the fact that a finite 
temperature is needed to melt the lattice in three dimensions. The melting of the flux-line lattice 
has been intensively discussed in connection with the irreversibility line in the phase diagram of 



high temperature superconductors [12| [13|, |J, 15], [§]. The question addressed by this paper is 



what happens if the random background potential, static in the first approximation, is allowed to 
fluctuate slowly in time due to the thermal hopping of pinning centers. We shall study this problem 
for a three-dimensional triangular lattice with a dissipative dynamics, weakly pinned by hopping 
point defects. 

The difference between static and thermal disorder in three dimensions stems from the fact 
that the statistical average of the static random pinning force over a large distance grows faster 
with the distance than the statistical effect of the random force due to thermal fluctuations. The 
hopping of pinning centers softens, through time average, this devastating large-distance effect 
of the static disorder. At high hopping rate its effect on the elastic lattice is equivalent to the 
thermal noise. We will argue that the translational order in the lattice should be restored when 
the rate of hopping attains some critical value and will discuss possible implications of that for 
superconductors. This phenomenon could also occur in charge density waves |TB| and Wigner 
crystals in semiconductors, atomic monolayers on imperfect surfaces [0, and other lattice 
systems where hopping of pinning centers may exist. 

In a continuous approach, deformations of the lattice are described by the displacement field 



u(r). The free energy of the system is [T9| 



~ / d 6 r [(C n - C m ){d a u a f + C m (d a upf + C u (d z u a )' 2 } - I d :i ru a f a 



2. 

where a, (3 can be x or y\ Cn, C 44 , and C 66 are the elastic moduli of the triangular lattice; f (r, t) is 
the force associated with random pinning. It should be noted that the replacement of the random 
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potential by the random force term in Eq.(l) is only true for small displacements as it fails to 
account properly for the periodicity of the lattice and the potential nature of the interaction with 
the random background || ||, [7|, j|, ||, p7D[| . For our purpose, however, the random force model should 
be sufficient because it gives the correct estimate of the correlation length. The static Gaussian 
noise is equivalent to < / a (q)/ J s(q') >= (2ir) 3 W6 a /35(q + q') for the Fourier transform of f(r). 
Here W characterizes the strength of pinning. The hopping of pinning centers can be introduced 
by assuming that < / a (q, t)fp{q[, t') > is proportional to exp(— r|£ — t'|), where T is the rate of 
hopping. Since pinning sites are independent, both W and T should be independent of q [pOfl . This 
gives 
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T 2 + u 2 

At not very low temperatures the dynamics of the flux-line lattice is dissipative [|J, 

du 5F 
^ dt Su 

with F given by Eq.(l) and r\ being the drag coefficient. Substituting into this equation 



u(r,t) 



one obtains for the Fourier transform of u 



(2vr) 4 U ^ q '^ e 



(2) 



(3) 



(4) 



iwqu a + (Cn - C 66 )(q-u)g Q + (C m q]_ + C u q 2 z )u a = f a (q,uj) 



(5) 



where the elastic moduli, in general, depend on q; q\ — q 2 + q 2 . The solution of this equation is 



w Q (q, oj) = (iuji] + Cuq]_ + C 44 ^) -1 — - + {iujrj + C m q{ + C i4 q 2 z )~ l 



fa 



(6) 



2 < [u(r,t)-u(0,t)] 2 > 



(7) 



To determine the correlation length we shall compute 

B<r) = i 

where a is the lattice parameter. This can be accomplished by substituting here Eq.(4) with 
u(q, u) given by Eq.(6), and averaging over random forces with the help of Eq.(2). Working out 
delta-functions one obtains 

1 1 



B(r) = I - , . — -[1 - cos(q-r) 



2a 2 J (2vr) 4 T 2 + lu< 



+ 



u;V + (C liq l + C m 2 z ) 2 u 2 r, 2 + (C m q\ + C m 
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Further integration over oj gives 

B(t) = ^/^[l-cos(q-r)][(rr / + C'ngi + C 44 g 2 2 )- 1 (C'ngi + C 44 g 2 2 )- 1 

+ (Trj + Cwqi + Cuqlr^Cwqi + Cuq!)- 1 } . (9) 

The effect of hopping now becomes apparent. It softens the small q singularity under the integral. 
When T goes to zero, Eq.(9) reduces to the well-known expression for B in the case of a static 
disorder [p], |J. 

We shall now perform the final integration over q in Eq.(9). At small V the integral is dominated 
by small q. To shorten the formulas we shall assume that C\\ » Cqq, which is true for flux lattices 



at small q | 21| . Introducing 



D (2,2, °66 2 \ 

R = [ x +y + cZ z ) 

R = l^a 2 C\i 2 C% 2 /W 

* - (Wf (10) 

and rescaling variables under the integral, one obtains after a series of standard integrations 

^^-tH-M)]} • <»> 

Here Rq is the Larkin-Ovchinnikov correlation length Jl| due to the static disorder, while Rh is the 
length of the diffusion of the unpinned flux-line lattice during the time 

Formally, the translational correlation function given by Eq.(ll) diverges at R—>oo only when 
Rh^-oo, that is, at T^O. In that limit Eq.(ll) coincides with the Larkin's result, B = R/Rq. At 
any non-zero T, B approaches a constant 2Rh/Ro at R^oo. This result should be taken with 
caution, however. Indeed, Eq.(ll) does not provide the correct answer for translational correlations 
at large distances as it has been derived in the random force approximation, without taking into 
account the periodicity of the lattice. The answer for the random potential problem is not known. 
It is believed, however, that the random force approximation provides a good estimate of the 
translational correlation length, R t . The latter is the distance at which the sites of the deformed 
lattice can no longer sustain one-to-one correspondence with the sites of the perfect lattice. It can 
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be obtained as a solution of the equation B(R t ) = 1. At T^O, B = R/Rq, and, thus, R t = Rq. 
As T increases, Rt = Rq remains a good approximation at Rh » Ro- The significant change in R t 
occurs when Rh becomes comparable with Rq. The translational correlation length grows rapidly 
as Rh decreases below R and becomes infinite at Rh = -Ro/2- One should expect, therefore, that 
the long range translational order becomes restored at 

r = (12) 

The limit of Rh large compared to Rq corresponds to small T and/or small 77 [fZ^]. In that 
limit the lattice adjusts to the evolution of the pinning potential fast enough to make this situation 
similar to the case of the static disorder. In the opposite limit of Rh small compared to Ro, the 
pinning potential changes fast enough to make its effect on the lattice similar to the effect of thermal 
fluctuations. The latter effect is equivalent, through the Einstein relation, to the random force with 
the correlation function 

< / a (q, u)f p (q', lo') >= (2n) 4 2Tr ] 5 a/3 6(q, + q')5(uj + «/) . (13) 

Repeating all steps leading to Eq.(9), one obtains 



1 1 

+ 



(14) 



C n ql + C u q 2 z C m q]_ + C u q 2 _ 

for the thermal disorder. This expression is finite at all R and allows to estimate the melting 
temperature through the Lindemann criterion |fl2| , |i~3| , 14, 15, B(a) = c|, where cl < 1 is the 



Lindemann number. Equations (9) and (14) become identical in the limit of large V (V » 
CggQbz/v) ^ one replaces W/T in Eq.(9) by 77T, which represents the transition from hopping to 
thermal noise. Here qsz = (87r) 1 / 2 /3 1//4 a is the radius of the Brillouin zone. It is easy to see that as 
long as Rq remains large compared to a, that is, in the weak pinning regime, T c of Eq.(12) is always 
small compared to the hopping rate, T ~ C^q^z/rj, at which hopping crosses to the thermal noise. 

For a superconductor Cqq = B c2 B^ /327rA 2 |l|, |J, where £ and A are the coherence length and 
the penetration length respectively. The dissipation coefficient r\ equals [^3| B c2 §q/ p n c 2 a 2 , where 
$0 is the flux quantum, p n = 47rA 2 z/ e /c 2 is the normal state resistivity (y e being the normal electron 
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collision frequency), and a = (2 1/,2 /3 1//4 )(<3>o/-B) 1//2 is the flux lattice spacing. Substituting all that 
into Eq.(12), one obtains a remarkably simple expression for T c : 

Let us now study the question whether the corresponding rate of hopping may exist in super- 



conductors. Measurements of translational correlations in YBCO and BSCCO [M, |25| suggest that 
at a low field the ratio £/-Ro may be as small as 10~ 3 . With v e ~ 10 13 s _1 , it requires T as high as 
10 7 s _1 . In principle, hopping of defects may be present at zero temperature due to quantum tunnel- 
ing, but the rate should be small. Thermal hopping of atomic defects in high-T c superconductors, 
at a sufficiently high temperature, can be much stronger. Assuming that V follows the Arrhenius 
law, T = z/oexp(— U/T), with v$ of order of the Debye frequency, 10 13 s _1 , we find that at 77 K the 
barrier U that would provide the critical hopping rate must be about 1000K. Barriers as low as 870 
K, presumably related to hopping of isolated oxygen atoms in CuO planes, do exist in YBCO, see, 



e.g., Ref. pra and references therein. In copper oxygen superconductors without twin boundaries, 



BiSCCO in particular, oxygen vacancies can be the major source of pinning [p7 |. The significance of 
their hopping could be established, e.g., through measurements of the magnetic relaxation induced 
by a sudden change in the pressure of oxygen. If a significant relaxation is detected, it would be 
quite conceivable that the phenomenon suggested in this paper could take place in the low-field 
high-temperature (large T) part of the phase diagram. The restoration of the long range order in 
the flux-line lattice due to the hopping of pinning centers would be effectively equivalent to the 
depinning of the lattice. It would result in the critical current, j c , dropping to zero because j c goes 
down as a some power of Rt 0. Some additional experimental effort may be needed, however, to 
distinguish between the effect of hopping and the conventional thermal depinning or melting of the 
flux-line lattice. 

This work has been supported by the Department of Energy under Grant No. DE-FG02- 
93ER45487. 
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